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Abstract. A binary mixture made of heavy and light inelastic hard spheres in the free cooling state is considered. First, the 
regions in the parameter space where the partial granular temperature of the heavy species is larger (or smaller) than that of 
the light species are analyzed. Next, the asymptotic behavior of the mean square velocity ratio in the disparate-mass limit is 
investigated, assuming different scaling laws for the parameters of the mixture. It is seen that two general classes of states are 
possible: a "normal" state and an "ordered" state, the latter representing a strong breakdown of energy equipartition. 
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As is well known, granular fluids are intrinsically out of equilibrium. This clearly manifests itself in the case of 
granular mixtures, where the equipartition of energy is broken down (even in homogeneous and isotropic states), as 
shown by theory 001, simulations [3], and experiments |4]. In this work we consider a binary mixture made of heavy 
(h) and light (£) inelastic hard spheres in the free cooling state. By using kinetic theory tools, one can derive a tenth- 
degree algebraic equation 1 1] whose solution gives the mean square velocity ratio = (v?)/(v?) = 7/,m//7)m/, (which 
plays the role of an "order" parameter) as a function of the control parameters of the problem (the densities, the three 
coefficients of normal restitution, and the mass and size ratios). We first investigate the regions in the parameter space 
where the partial granular temperature of the heavy species is larger (Tf, > Tg) or smaller (Tf, < T() than that of the light 
species. We next analyze the asymptotic behavior of <f> in the disparate-mass limit (me/m/, — > 0), assuming different 
scaling laws for the coefficients of restitution, size ratio, and concentrations. We observe that it is possible to distinguish 
two general classes of states: a normal state (where 0^0) and an ordered state (where ^ 0), the latter representing 
a strong breakdown of energy equipartition in which the temperature ratio Tf,/T( diverges. These two classes can be 
further resolved. Thus, the normal state can be partitioned (T^/Tf — finite), mono-energetic h (T/,/T( — > °o), or mono- 
energetic I (Th/T( — > 0). Moreover, if the concentration of heavy particles vanishes, the ordered state can become 
super-ordered (0 — * °°). The phase diagrams corresponding to this rich phenomenology are presented. In particular, in 
the limit of vanishing concentration of heavy particles, the partitioned phase is separated from the ordered and from 
the super-ordered phases by two critical lines. This extends a previous analysis made in the case of an impurity particle 
immersed in an inelastic gas 



Let us consider the homogeneous cooling state of a binary mixture of two granular species (h and £) characterized by 
mole fractions x/, = «/,/« andxf = 1 — jc/„ masses andmf < m/„ diameters a%, <J(, and Ohi = (o), + o>) /2, coefficients 
of normal restitution a/,/,, c%, and 0£/^, and pair correlation contact values Xhh> Xu< an d Xhl- This gives 10 independent 
parameters. However, in appropriate reduced units the number of control parameters reduces considerably. The set of 
two Boltzmann equations are 



INTRODUCTION 



BASIC EQUATIONS AND CONDITIONS FOR T h > T e 



Multiplying by v 



dtfh{v) = Jhh[v\fhJh] +Jhl[v\fh,ft], 
• 2 and integrating over velocity one gets 
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d,{v 2 h ) = + d t {vj) = -(C« + C«)(vl>, £ 
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The diagonal terms and Q( are positive definite and represent the inelastic cooling rates of species h and £, 
respectively. On the other hand, the cross term §^ represents the "thermalization" rate of species h due to collisions 
with particles of species I and can be either positive or negative, depending on the state of the mixture. Analogously, 
the term Qt, represents the thermalization rate of species £ due to collisions with particles of species h. 

The four rates £y are complicated nonlinear functionals of the unknown distribution functions //, and fg. However, 
reasonably good estimates can be obtained by using the Maxwellian approximation 

fi{\) = m (m ( /2^7i) 3/2 exp (-m ( v 2 /27-) , 7- = m,(vf)/3. (3) 

This approximation yields £y = CO^y, where (0 = 4^«^a^(v^)(l + a^)m^/3(m/, + m.() is an effective collision 
frequency and the dimensionless cooling and thermalization rates are [ 1] 



Zhh(<i>)=XhV<j>Ph, = x e~T~7~r~ ( 1 - M + 00 - x ) i ( 4 ) 

1+00 V P. 
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Here, = (v 2 ,) / (v^) is the mean square velocity ratio, fx = mi/im^ + ntg) is the mass ratio, 0o = (1 — /(l + OE^) 
is a measure of the cross coefficient of restitution, and the coefficients 

g _ 1 +<t>Q 1 ~ a hh Xhh ( Qh V g l + <hI-"u Xu ( °e\ 2 (6) 



essentially measure the self coefficients of restitution and are inversely proportional to ji. The evolution equations (|2ji 
yield (0~ l dt(j) — (^a + ^th — %hh — <^/if)0i so that the homogeneous cooling state of the mixture is characterized by the 
stationarity condition 

&*W+§M(«)-§«W-&kW=0. (7) 

This equation can be transformed into a tenth-degree algebraic equation for 0. Its physical solution determines 
as a function of the five control parameters of the problem, namely {xf,, (Xm, OCgg, <XM,mi/m$ or, equivalently, 
{x h ,p h ,p e ,(j)o,fl}. 

Equipartition of energy implies that = ju/(l — fl). Of course, this happens when all the collisions are elastic (i.e., 
j3/, = fie = 0o = 0)- When collisions are inelastic, the most common situation is > — /x), i.e., the heavy particles 
have a larger (partial) temperature than the light particles, 7), > Tg |3j|4j]. However, this is not always necessarily the 
case. To clarify this point, it is convenient to define the parameters Kj = lij3,(l + 0o)/0o> i = which essentially 
measure the degree of inelasticity of the i-i collisions, relative to the inelasticity of the h-l collisions. A careful analysis 
of Eq. |7} shows that, for a given value ji < j, one has Tf, > Tg when one of the following three alternative cases 
occurs: (i) K h > 7T*(ii), Kg > X|0), and x h < x* h (K h ,Kg,n); (ii) K h < £|0), Kg < Kj[(ji), and x h > x* h (K h ,Kg,^); 
(iii) K h < K* h (n) and Kg > K^(jx). Here, 



The left panel of Fig.^sketches the regions in the Kj-Kn plane where 7), > Tg or 7), < 7j for a given mass ratio. The 
right panel shows the temperature ratio versus concentration for a few representative cases. 



THE DISPARATE-MASS LIMIT PHASE DIAGRAMS 

We are now interested in the limit mg/mf, — ► => ji — > 0. The question we want to address is, how does the solution 
of Eq. Q behave in the limit ji — > 0? There are two distinct possibilities: either lim^o = or lim^o 7^ 0. Thus, 
can be viewed as an "order" parameter. The case lim^o 0=0 means that the heavy particles move much more 
slowly than the light particles; we will say that the mixture is then in a normal state. The alternative case lim^o 7^ 
implies that the speed of the heavy particles is typically comparable (or even larger) than that of the light species; we 
will call ordered to such a state. 
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FIGURE 1. Left panel: Plane Kf, vs Kg split into four quadrants by the lines AT/, = K'^(fl) and = for = ^; in the 

top-right quadrant, 7), > 7f if x/, < ^{K^K^fl), as represented by the inset, where x* h is plotted as a function of K( for /l = ^ 
and = 2; in the bottom-right quadrant, 7), > 7^ for any concentration; in the bottom-left quadrant, 7), > 7) if x/, > £{Kh,Ki,ll), 

as represented by the inset, where x^ is plotted as a function of Ki for /x = J and iC), = 0.5; finally, in the top-left quadrant, 
7), < 7^ for any concentration. Right panel: Plot of the temperature ratio 7/,/7f versus x/, for (Xm = 0.8, m/, = 2m^, o>j = Ci, and (a) 
(Otok.ote) = (0.6,0.6), (b) {a hh ,a u ) = (0.9,0.6), (c) (a hh ,a u ) = (0.9,0.9), and(d) (a hh ,a u ) = (0.6,0.9). 



In order to make a finer classification, let us introduce a "critical exponent" 77 such that <p ~ fi^ as — * 0, so 
77 > corresponds to normal states and 17 < correspond to ordered states. If T7 = 1, the equipartition of energy is 
weakly broken in the sense that 7/,/Tf = finite, as in the elastic case. This corresponds to a subclass of the normal 
state that we will call partitioned. An example of this situation occurs when the self-collisions are elastic and the 
cross collisions are quasi-elastic, with 1 — ay, ~ fl- On the other hand, if 77 > 1, then (j) decays even more rapidly 
than in the elastic limit. Again, this is a normal state (<f> — > 0) but is not partitioned (7)j/7> — > 0) since all the energy 
is carried by the light particles. This case will be referred to as a mono-energetic I state, an example of which taking 
place when only the h-h collisions are inelastic. Still in the normal class, the case < f] < 1 represents a breakdown of 
equipartition intermediate between that of the partitioned state and that of the ordered state: while the heavy particles 
move much more slowly than the light ones, they carry all the energy of the mixture. We will refer to this situation 
as a mono-energetic h state. A simple example corresponds to mixtures of disparate sizes when only the £-£ collisions 
are inelastic. Thus, three subclasses (mono-energetic I, partitioned, and mono-energetic h) can be distinguished within 
the class of normal states. As for ordered states, the typical case corresponds to tj = 0, i.e., <p = finite in the disparate- 
mass limit. This happens, for instance, when only the cross collisions are inelastic, regardless of the values of the mole 
fractions. In particular, the state is ordered both in the "Brownian" limit xj, — » and in the "Lorentz" limit X( — > 0. 
In the former case, 0^0 implies that the dynamics is not Brownian since the solute Brownian particles move with a 
speed comparable to that of the fluid particles. In the latter case, the system does not behave as a Lorentz gas since 
the heavy "scatterers" are never at rest. An even stronger breakdown of energy equipartition occurs if rj < 0, implying 
not only ^ (ordered state), but (j) — > 00. We will say that in this case the system reaches a super-ordered state. A 
necessary condition for the existence of this situation is the Brownian limit Xf, — > 0. A super-ordered state appears, for 
instance, if only the £-£ collisions are inelastic and Xf, ~ fl, so that the total masses of both species are comparable. 



TABLE 1. Possible classes of states in a binary granular mixture in the disparate-mass limit 
ntf/ntf, ~ /x — > 0. The mean square ratio <p = (v?)/(v?) scales as <j> ~ /i 7 '. 
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Table^describes the five possible subclasses (or phases) and the corresponding representative examples. In general, 
depending on the regime of values for the control parameters a,- ; and Xh, the asymptotic homogeneous cooling state of 
a disparate-mass binary mixture belongs in one of those subclasses. In order to investigate the regions in the parameter 



space corresponding to each phase, let us assume scaling laws of the form 

p h ~H- 1+a », Pt~n- 1+a *, *h~H Cl \ (9) 

where the exponents a%, ag, b, and c/, are non-negative. The exponents a/, and b measure, respectively, the inelasticity 
of the h-h and h-t collisions (i.e., 1 — ~ jJ," 1 ', 1 — ~ ju*), while measures the inelasticity of the £-£ collisions, 
as well as the possible dependence of the size ratio on the mass ratio: (1 — au)(<Je/ 'Oft) 2 ~ fl at . If any of the three 
types of collision is elastic, the corresponding exponent takes an infinite value. Conversely, if the collision is inelastic 
then the exponent vanishes. A finite value of a/, or b implies that the associated collision is quasi-elastic. In the case 
ag — finite, either the £-£ collisions are quasi-elastic or the size ratio scales with the mass ratio, or both. For finite 
concentrations of the h species one has c/, = 0, while c/ 7 > in the Brownian limit (x/j — > 0). From Eqs. and (0, the 
scaling behaviors of £y become 

&ft ~ jll-W^^ ^, „ maxlM^M^^maxlM , \^ |}, (10) 



fe~M" 1+a ^ ^^M c "max{M ,M 77/2 }max{M ,M" 1+ MM" 1+7? |}- (H) 

In the above equations, | • • • | indicates that the corresponding term has a negative prefactor. 

Let us first assume that the mole fraction Xf, is kept finite in the limit jj. — > 0, so that c/, = 0. In that case, condition 
|7} is fulfilled in the limit jj. — > when any of the five conditions described in Table [2] is satisfied. The associated 
constraints and expressions for the exponent rj are also included in Table [2] It is easy to check that in region B one 
always has rj > 1, while in region D one always has < TJ < 1. In region E, TJ = 1 if Z> > 1, while < tj < 1 if 
< b < 1. As for region A, < 77 < 1 if an < a/, + j, while 77 > 1 otherwise, this latter case being only possible when 
b > i. Likewise, in region C one has < 77 < 1 if at, > b — | and 77 > 1 otherwise, this latter case being only possible 
when b> j. The phase diagrams are shown in Fig.|2] 



TABLE 2. Possible cases in the disparate-mass limit at finite concentration (c/, = 0). 
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FIGURE 2. Regions in the plane ci( vs a/, corresponding to each phase. The thick lines separating the phases mono-energetic 
£ and mono-energetic h correspond to partitioned states. The thick line at ag = represents ordered states. The whole plane is 
occupied by ordered states if b = (not shown). 

Now we assume that the concentration of h particles behaves as Xf, ~ fl c '' in the limit jj. — > 0. The total mass of the h 
species is larger than that of the £ species (i.e., x^nif, ^> Xgm£) if c/, < 1, while the opposite happens if c/, > 1. This limit 
adds a new parameter cj, > to the parameter space, so that the phase diagrams in Fig.|2]are modified. Let us restrict 
ourselves to b — 0, i.e., the cross collisions are kept inelastic in the limit jtl — > 0. As indicated in the caption of Fig. 
|2j if b = the state of the mixture is always ordered (77 = 0) in the case of finite concentrations (q, = 0). However, 
this uniform situation changes dramatically if cj, > 0. With b = 0, it is easy to realize from Eqs. (II Oi and (lilt that 



is never negligible versus ^/,/,. That means that the value of the parameter ah (i.e., the degree of inelasticity of the 
h-h collisions) will not play any role, so the parameter space reduces to ai and c/,. The possible cases are described 
in Table|3] Cases Gi and G2 correspond to maxf/i ,/^/ 2 } = jit* 1 ' 2 and max{^°,^ T '/ 2 } = jU°, respectively. While Gi 
defines a whole region in the parameter space with a well-defined (negative) expression for the exponent 77, G2 defines 
a line along which the value of 77 is not unique. A similar situation occurs with case H. The phase diagram is shown 
in Fig. [5] We observe that 77 changes discontinuously from 77 = to 77 = 1 when going from the ordered region to the 
partitioned region through line H. Even more discontinuous is the transition from the super-ordered region (77 < 0) 
to the partitioned region (77 = 1) through line G2. These are critical lines along which the value of 77 is not uniquely 
determined by the values of ai and c/,. Let us analyze these critical lines with some detail. 



TABLE 3. Possible cases in the disparate-mass limit with Xf, — » and 
a ht < 1 (b = 0). 
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FIGURE 3. Regions in the plane ci( vs c/, corresponding to the phases present in the limit of vanishing concentration of heavy 
particles when 1 — a^i is kept finite (b = 0). The thick line separating the regions super-ordered and partitioned (line G2) and the 
thick line separating the regions ordered and partitioned (line H) are critical lines on which the state can be partitioned, mono- 
energetic h, or ordered, depending on the value of the parameters ft, (line G 2 ) or y (line H). 



Along the critical line G2 the leading terms are ^ and so that condition Q becomes 



The solution to this equation in the limit ji — * is 



J" 



0(l + 0o) 



(12) 



1 



rife if Z^ 1 ' 
J3 2 -1 ifft>l. 



(13) 



Therefore, j3<- plays the role of a control parameter and fa = 1 is its critical value. If fa < 1 the state is partitioned, 
if fa = 1 it is mono-energetic h, and if fa> > 1 it is ordered. This situation has been analyzed in detail in Ref. @]. 
Since Cf, > 1, the total mass of the h species is negligible versus that of the I species. Moreover, the £ particles are not 
affected by the presence of the h particles. 

Along the critical line H, Xf, ~ fi (i.e., c/, = 1), what means that both species have comparable total masses. In 
addition, the l-i collisions are not too inelastic (or the size ratio O^/cT/, is sufficiently small), in the sense that 



(1 — a^)((7^/<7/j) 2 ~ with a/ > 1. In that situation, the leading terms are t,^ and Therefore, condition Q 
yields 

1 

whose solution is 



Here, y= (x/,//x)0o/(l + 0o) = (x/,/jti)(l — a/,f)/2, which represents the control parameter in this case. 

CONCLUSIONS 

In this paper we have shown that in the homogeneous cooling state of a granular binary mixture the mean square 
velocity ratio = (v?)/(v|) and the temperature ratio Tf,/Tf can take widely different values depending on the 
parameters of the system (coefficients of restitution, concentrations, and mass and size ratios). Typically, 7), > Ti 
if (i) the self-collisions {h-h and £-£) are sufficiently more inelastic than the cross collisions (h-£) and the concentration 
of the h particles is sufficiently low; (ii) the h-h and £-£ collisions are sufficiently less inelastic than the h-£ collisions 
and Xf, is large enough; (iii) the h-h collisions are not too inelastic but the £-£ collisions are sufficiently inelastic. 

In the disparate-mass limit fi — > 0, the breakdown of equipartition can become extreme with divergent, finite, or 
vanishing limits of the mean square velocity and/or temperature ratios, depending on the scaling behavior of the 
parameters. This gives rise to the five classes of states (or phases) described in Table [0 ranging from the mono- 
energetic £ state (JhjTt — ► 0) to the super-ordered state (0 — * °°). The associated phase diagrams are presented in Figs. 
[2]and[3] If the cross collisions are inelastic (a^ < 1), the state is always ordered (0 = finite). As a consequence, in 
this case there is neither Brownian dynamics (when X}, — > 0) nor Lorentz dynamics (when X( — > 0). A partitioned state 
{Ti,/T( = finite) is only possible if the three types of collisions are sufficiently quasi-elastic, while a super-ordered state 
(0 — > oo) is only possible in the Brownian limit (xj, — + 0). In that limit, there is no mono-energetic £ state and there 
exist critical lines in the phase diagram where the state can be partitioned, ordered, or mono-energetic h. 

The analysis performed in this work has been restricted to the homogeneous cooling state, but it can be easily ex- 
tended to mixtures heated with a white-noise forcing 1 2 ] . In that case, condition |7} is replaced by [§/ja(0) + — 
<^«(0) — E,eh(<P) = 0. As a consequence, the phase monoenergetic £ is suppressed. 

ACKNOWLEDGMENTS 

A.S. acknowledges partial support from the Ministerio de Educacion y Ciencia (Spain) through Grant No. FIS2004- 
01399. 



REFERENCES 

1. V. Garzo, and J. W. Dufty, Phys. Rev. E 60, 5706-5713 (1999). 

2. A. Barrat and E. Trizac, Granular Matter 4, 57-63 (2002). 

3. M. Alam and S. Luding, Granular Matter 4, 139-142 (2002); S. R. Dahl, C. M. Hrenya, V. Garzo, and J. W. Dufty, Phys. Rev. 
E 66, 041301-1-10 (2002); A. Barrat and E. Trizac, ibid. 66, 051303-1-13 (2002). 

4. R. D. Wildman and D. J. Parker, Phys. Rev. Lett. 88, 064301-1^1 (2002); K. Feitosa and N. Menon, ibid. 88, 198301-1-4 
(2002). 

5. A. Santos and J. W. Dufty, Phys. Rev. Lett. 86, 4823-4826 (2001); Phys. Rev. E 64, 051305-1-16 (2001). 



Xf, 0() - 



0(1 + 00) JU 1 + 0O : 



if 7 < 1 , 
if 7=1, 



l -y 

1+00 

ffi(y-l) if y> 1 



(14) 



(15) 



